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Abstract 

With the discovery of a modest size for the mixing angle ^13 ^ 9° by the Daya Bay collaboration at >5 

fvj sigma ([!]) the science of neutrino oscillations has shifted to explicit demonstration of CP violation and 

^-H precision determination of the CP phase in the 3-flavor framework. Any additional contributions from new 

i^ physics to the oscillation channel z/^ — > Ve could be uncovered by multiple constraints in the (0i3,(5cp) 

. parameter space. In long-baseline experiments such constraints will require examination of the oscillation 

^2 strength at higher L/E where the effects of CP violation will be large. For the fixed baseline of 1300 

^5 km for the Long-Baseline Neutrino Experiment (LBNE, Fermilab to Homestake), it will be important to 

Q\ examine oscillations at low energies (< 1.5 GeV) with good statistics, low backgrounds, and excellent energy 

'"^ resolution. The accelerator upgrades in the Project-X era have the potential to offer the beams of the needed 

I— I intensity and quality for this advanced science program. In this paper we examine the event rates for high 

^ intensity, low energy running of Project-X and the Fermilab Main Injector complex, and the precision in 

2^ the (013, 5cp) space. In this paper we have examined the baseline distance of 1300 km in detail, however we 

1^ point out that much longer distances such as 2500 km should also be exmained with a beam from FNAL in 

r^ I light of the new understanding of the neutrino mixing. 

We find that the best way to obtain multiple constraints in the neutrino sector is to perform high statistics 

►^ experiments with low energy neutrino beams over long distances. For oscillation physics at low energies the 

^^ charged current cross section is dominated by quasielastic scattering. For the quasielastic final state any 

^iJ large detector capable of measuring single lepton final states is adequate. In the following we have used the 

'sj" water Cherenkov detector since it can be built to have the target mass (^ 200kTon) needed to obtain the 

(T^ required statistical precision and, as a result of the work for LBNE, the performance of the detector at low 

^< energies is well established. 



X 



I. INTRODUCTION 

The Daya Bay Reactor Neutrino Experiment has measured a non-zero value for the neutrino 
mixing angle ^13 with a significance of 5.2 standard deviations [Ij. The value determined by using 
a rate-only analysis of the first data set from the experiment is sin^ 2^13 = 0.092 it 0.016 it 0.005. 
The error currently is dominated by statistics. 

This determination of a modest size for the third mixing angle (^13 ~ 9°) in the neutrino sector 
allows us to plan for the future program of precision neutrino physics. With this value for the third 
mixing angle, a rich new system has become available for the physics of neutrino oscillations, and 
leptonic CP violation. 

Now that the parameters of neutrino mixing have been measured it will be necessary to exploit 
the new system to find indications of departures from the Standard Model with massive neutrinos 
(z^SM). These departures could come from additional low energy interactions or unexpected mixing. 
It is not possible to predict the exact nature of new physics or investigate every model for new 
neutrino physics that is now being discussed in the theoretical literature. Nevertheless, it is possible 
to sketch some generic attributes regarding the design of a new experiment that is likely to yield 
indications of new physics. These necessary attributes would be based on physics that is known, 
physics that has already been excluded and reasonable extrapolations of scientific strategies that 
have worked in the past. 

An almost complete history of neutrino oscillation experiments (except for the Daya Bay ex- 
periment) can be found in figure [11 There are some important lessons to be learned in this figure. 
Experiments designed for high energy operation probe mixing in the region of Am^ > leV'^. Many 
such experiments were carried out with sensitivity to smaller and smaller mixing angles with limited 
success. Eventually mixing was found to be large and the correct mass range to be much smaller 
than 1 eV'^. Previous results as well as the direct limits on neutrino mass from tritium beta decay 
(m,jy^ < 2 eV at 95% C.L.) and astrophysics {J2i ^^Vi <0.4 - 1 eV) have put severe limits on physics 
at L/E < 1 km/GeV [2j. Furthermore, any remaining issues such as the LSND/MiniBoone/reactor 
anomalies are best resolved by short baseline accelerator experiments. As we examine the oscil- 
lation phenomena in greater detail we find that large interesting effects are always at high L/E. 
Moreover, most of the precision measurements in neutrino physics have been performed at low 
energies. Experiments that have yielded the best results are at high L/E > 500 km/GeV, or for 
any fixed distance, the lowest energy neutrinos have provided the best science. 

In the following, we first calculate the oscillation probability for f^ — )• fe as a function of energy 
using the best values for the oscillation parameters. We will then calculate the energy spectrum 
of neutrinos that could be produced by various beams from Fermilab. Finally we calculate the 
precision in the (^13, Scp) space for an example running condition, and comment on constraints on 
other parameters that could result from the data. 

For oscillation physics at low energies the charged current cross section is dominated by 
quasielastic scattering as shown in figure [2| For the quasielastic final state any large detector 
capable of measuring single lepton final states is adequate. The water Cherenkov detector can 
achieve the target mass needed to obtain the statistical precision at low energies with high effi- 
ciency. In the following we will therefore assume a 200 kTon fiducial water Cherenkov detector for 
the target mass. Furthermore, as a result of the work for LBNE the performance of the detector 
at low energies is well established. 



II. Vf, ->■ Ve OSCILLATIONS 

The oscillation mode with the richest scientific content is the appearance mode t'^ — t- z^e (or 
t'e — ^ t^fi) and its antineutrino counterpart. The oscillation strength for this mode as a function of 
distance and energy contains contributions from all parameters of the neutrino mixing matrix and 
all three neutrino mass eigenstates. 

Assuming a constant matter density, which is a good approximation for the baselines in question, 
the oscillation of f^ — )• Vg in the Earth for 3-generation mixing is described approximately by the 
following equation [5] 
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FIG. 1: Compilation of neutrino oscillation limits and results. The filled areas are positive indications and 
lines correspond to 90% confidence level limits unless otherwise indicated. [3] 
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FIG. 2: Charged current neutrino total cross section as a function of energy for neutrinos (left) and antineu- 
trinos (right). The cross section of quasielastic interactions is also displayed. The error bars are the data 
and the lines are the current best model. It should be noted that below 1.5 GeV the total cross section is 
dominated by quasielastic interactions. [3] 
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where a = Am^i/ Amli, A = Am^^L/AE, A = 2VE/Am'^i, V = \/2Gpn(.. rie is the density of 
electrons in the Earth. Recah that Am^i = Am'^2 + A'tiIi- Also notice that AA = LGpne/v^ is 
sensitive to the sign of Atti^i. For antineutrinos, the second term in Equation 111 has the opposite 
sign, and the matter potential also has the opposite sign. The second term is proportional to the 
following CP violating quantity. 

Jcp = sin 6*12 sin 023 sin ^13 cos 6*12 cos ^23 cos ^13 sin 5cp (2) 

Equation[T]is an expansion in powers of a. This approximate formula is useful for understanding 
important features of the appearance probability: 1) the first three terms in the equation control 
the matter induced enhancement for normal mass ordering (mi < ?tt-2 < fn^) or suppression for the 
reversed mass ordering [m^ < mi < m2) of the oscillation probability above 3 GeV; 2) the second 
and third terms control the sensitivity to CP in the ~ 1 GeV range; and 3) the last term controls 
the sensitivity to Am^i at low energies. The first term (last term) is also proportional sin^ ^23 
(cos^ ^23)) and therefore is sensitive to the issue of maximum mixing in ^23 = vr/4. As previously 
explained [6], measurement of the spectrum of oscillated fg events will allow us access to all of 
these parameters in a single experiment with good control of systematics. 

The oscillation probability for f^ — )• I'e for both neutrino and antineutrino modes and for normal 
and reversed mass ordering is plotted in figure [3] using close to the current best known parameters 
[7] and rounding the Daya Bay measurement to the value sin^ 2^13 = 0.1. It is easy to see that the 
oscillation probability will be in the 5-10 % range at the first oscillation maximum between 2-3 
GeV, but at lower energies both the size of the probability and the effect of the CP phase will be 
dramatic. The lower energy oscillation effect is also relatively less affected by the mass hierarchy. 
A measurement of these probabilities across the energy range will certainly result in precise new 
information about the mechanism of neutrino oscillations. 



III. EVENT RATES FOR LBNE FOR VARIOUS PROJECT-X BEAM CONDITIONS 

The relevant energy range of oscillations over 1300 km can be seen from figure 1 to be from 
0.2 GeV to 4 GeV. Any flux above 4 GeV is not going to contribute much to physics, and events 
below 0.2 GeV will most likely have poor resolution. Neutrinos in the 0.2 to 4 GeV energy range 
are produced by pions of energy ~ 0.5 — 10 GeV. We shall assume that the proton energy available 
from Project-X and related upgrades will be in the 8 to 120 GeV range. The production spectrum 
of neutrinos in the laboratory frame is affected by the proton energy, the focusing efficiency of 
the target/horn system, and finally the kinematics of pion decay. In the center of mass of proton 
collisions, pion production is distributed in an approximate normal distribution as a function of 
rapidity and increases only logarithmically as a function of y^ oc ^JEprotons- In the laboratory 
frame this distribution of pions is boosted by a factor 7cm ~ ^/E^jroton- For higher proton energies, 
there are higher number of pions produced at high rapidity in the center of mass, and they are 
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FIG. 3: Appearance probability for i/^ — >■ i/g as a function of energy at a distance of 1300 km. The top plots 
are for neutrinos and bottom plots are for antincutrinos. The left side plots are for normal mass ordering and 
right hand side are for reversed mass ordering. The parameters used for these plots are Ato§2 — 0.00256^^, 
AmJ2 = 7.6 X 10-5eV'^, 6*23 = 7r/4, 6*12 = 34°, 6I13 = 9.2°. The blue curve in all cases is for Sep = and 
the red and green curves are for Sep — 7r/2 and — vr/2, respectively. 



boosted to higher energies in the laboratory frame. Both effects are such that in our proton and 
pion energy range, the production of pions at any given energy is roughly proportional to the 
energy of the protons and, obviously, the total current of protons, or the total proton beam power. 
If we assume that the target/horn system can be tuned to focus pions efficiently at any energy 
above 0.5 GeV, the neutrino yield is further affected by the kinematics of pion decay. The flux of 
neutrinos from pion decay in the forward region is proportional to 7^ oc E'^. The two kinematic 
effects - the production of pions as a function of proton energy and the pion decay kinematics ~ 
are such that the only way to boost the yield of neutrinos at low energies is with high proton beam 
power at low energies. 

In table |l] we have made a list of beam conditions that could be possible from Project-X and the 
Project-X upgrade at 8 GeV. Figure [4] shows the beam power available from the Main Injector as a 
function of energy. With Project-X the beam power from the Main Injector can be maintained at 
or above 2 MW over the range 60 - 120 GeV. This is because the decrease in energy can be (mostly) 
compensated by increasing the repetition rate. This trend continues as the energy decreases but at 
some point it is limited by the number of protons coming from the LIN AC. The power achievable 
at 30 GeV would be ~1.3 MW for the Project-X Reference Design. The Main Injector requires 270 
kW of incident 8 GeV beam power at 8 GeV to produce ~ 2 MW at 60 GeV. With the additional 
upgrade to the 8 GeV pulsed LIN AC, the 8 GeV power level could be increased to ~ 4 MW. In such 
a scenario, the Fermilab accelerator complex could produce multi-MW power at both 60 GeV and 
8 GeV simultaneously. The duty factor for any Main Injector operation would continue to remain 
small in the single turn extraction mode, however the duty factor at 8 GeV will be ~ 5 — 10% 
unless a ring is deployed to compress the beam further. The poor duty factor is not problematic as 
long as the far detector is deployed at depth in order to reduce backgrounds associated with cosmic 
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TABLE I: Beam conditions and power possible during the Project-X phase. An accumulator ring at 8 GeV 
could be used to improve the duty factor. [8] 



rays. The duty factor might be more important for the operation of the target/horn system. We 
will assume that a ring might be deployed to produce an appropriately short duty factor ~ 10~^. 

For our calculations of event rates and neutrino spectra, we will use operation at 60 GeV at 
2 MW and operation at 8 GeV at 3 MW. We have reduced the power assumption for 8 GeV 
running from 4 MW to 3 MW because some of the current (270 kW) will be needed for producing 
the simultaneous 60 GeV beam, and the rest might be needed for other experiments. We have 
calculated the beam spectra using a GEANT4 simulation of the LBNE beamline with magnetic 
horns with a current of 250 kAmps, 2 meter diameter decay tunnel with a length of 280 meters. 
The LBNE beamline was designed for high energy operation. There is currently no design to 
transport 8 GeV protons to the LBNE target. The 8 GeV beam will require either another beamline 
or substantial modifications to the current beamline design. Here we will not investigate these 
important technical issues regarding the beamline, but we will argue that simultaneous operation 
at 8 and 60 GeV is quite compelling and should be investigated. In the following we focus on the 
neutrino spectra and event rates. 

The muon neutrino and antineutrino spectra (without oscillations) for 8 GeV and 60 GeV beams 
are shown in figure [5] superimposed on the expected v^ — )• V(, oscillation probability. The event 
rate is calculated for the total muon neutrino (and antineutrino) cross section shown in figure 
[2j The two different energy spectra are shown to complement each other. The 8 GeV spectrum 
covers the low energy region where large CP phase effects exist while the 60 GeV spectrum covers 
the higher energy region where the matter effects will dominate. It should be remarked that the 
60 GeV beam also has similar numbers of events at low energies as the low energy beam, but 
the low energy beam is expected to have somewhat more rate and less backgrounds due to event 
mis-reconstruction. The beam contamination in these beams is shown in table HI] where the total 
event rate is tabulated for each component of the beam. The event rate after v^ disappearance is 
also shown for the muon neutrino component. A few comments are in order regarding this table: 

• The event rates have been calculated for the total cross section as in figure [2] In the 
next section we will use the tabulated water Cherenkov detector performance for extracting 
electron neutrino events and associated backgrounds. 

• The 60 GeV beam is very well tuned for the first oscillation maximum and consequently has 
a large effect due to muon neutrino disappearance. Almost 75% of the total muon neutrino 
events are calculated to disappear. This factor is smaller for the 8 GeV beam because of 
multiple oscillation nodes. 

• The neutrino contamination in the antineutrino beam is large ~ 26% for the 60 GeV beam, 
but it is much smaller for the 8 GeV beam. Nevertheless, the event rate for the antineutrino 
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FIG. 4: Proton beam power as a function of proton energy from the Fermilab Main Injector. Shown are 
current capabiUties labeled as NuMI. The recently funded upgrades (labeled as ANU) will increase the power 
to 550 kW at 60 GeV or 700 kW at 120 GeV. Project-X as currently conceived will allow beam power of 2 
MW at 60 GeV and 2.3 MW at 120 GeV. [S] 
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TABLE II: Total rate of events for the 8 GeV and 60 GeV neutrino and antineutrino beams. The 8 GeV 
(60 GeV) beam is assumed to have power of 3 MW (2MW). The running time is one year and the detector 
fiducial mass is 200 kTon at 1300 km from Fermilab. The event rate after the disappearance of muon 
neutrinos by oscillations is also shown for the muon neutrino and antineutrino components. 



running in the 8 GeV beam is much more suppressed compared to the neutrino running. This 
difference can be traced back to both the production rate of vr" and the neutrino/antineutrino 
cross sections at low energies. The large suppression of antineutrino rates at low energies 
is the additional reason for the complementarity between high and low energy running for 
LBNE. 
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FIG. 5: Spectra of event rates as a function of energy for 8 GeV (left) and 60 GeV (right) proton beams 
From Fermilab. The spectra are superimposed on the expected oscillation probability for normal hierarchy. 
Spectra are for the total charged current cross-section for muon neutrino (top) and antineutrinos (bottom). 
The beam is from Fermilab to Homcstake over a distance of 1300 km; the intensity for the 8 GeV beam is 
assumed to be 3 MW and for 60 GeV it is 2 MW. The detector size is 200 kTon fiducial mass. 

IV. ELECTRON NEUTRINO APPEARANCE AND PRECISION CP VIOLATION MEA- 
SUREMENTS 



We have tabulated the performance of the water Cherenkov detector using data analysis from 
Super-Kamiokande (SK). This performance has been verified by several independent checks on SK 
data. The performance of the 200 kTon LBNE water detector is expected to be similar or better 
because of higher pixelation and much better time resolution of the photo-multiplier tubes. Work 
continues to fully simulate and reconstruct events in the 200 kTon water Cherenkov detector. A 
recent eye-scan of the Monte Carlo events indicates considerable room for improvement beyond the 
Super-Kamiokande based performance by including low ring- multiplicity events [9j. Nevertheless, 
in the following we have used the SK based signal efficiency and background rejection since it is 
the most conservative estimate. 

The calculation is performed by first using the total charged current and neutral current event 
spectrum as a function of neutrino energy and using the efficiency for selecting events reconstructed 
as single electromagnetic showers with no signatures of additional particles including decays of 
muons that might be below Cherenkov threshold. This calculation includes the expected energy 
resolution and smearing of both charged current and neutral current events. After obtaining the 
spectrum of these events using tabulated efficiency factors, we use a tabulated event identification 
likelihood efficiency (LL) as a function of reconstructed energy. This likelihood performance was 
also tabulated using SuperK Monte Carlo that has been tuned to atmospheric neutrino data. 
The likelihood performance tabulation includes the efficiency for charged current electron neutrino 
signal events and both neutral current and charged current muon neutrino events. The likelihood 
performance can be tuned to obtain high signal efficiency (for example, from 80%LL to 40%LL) 
and different levels of background rejection [11]. 

We have used the parameterized performance for the water detector using the 60 GeV and 8 
GeV beams to calculate the reconstructed spectra with the proper resolution and efficiencies in 
figures [6] and [Tj For the 60 GeV running we have used a likelihood performance that has 40% 
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FIG. 6: Event spectra using a 200 kTon water Cherenkov detector at Homestakc for a neutrino beam 
from FNAL. We have assumed 2 MW of power at 60 GeV for 5 years of running for neutrinos (top plots) 
and antineutrino modes (bottom plots). The left plots are for normal mass hierarchy and right plots are 
for inverted mass hierarchy. We have put sin 20i3 = 0.1. The curve with error bars in each plot is for 
Sep = — 7r/2 and the red and blue curves are for Sep = 0, tt/2, respectively. The integral numbers of events 
are shown in the legends. The integral was computed up to the arrow in the figure. 



efficiency for signal in figure [6J The lower efficiency cut reduces the neutral current background 
above 2 GeV to be approximately the same as the beam contamination of electron neutrinos and 
provides a spectrum with a large signal to background ratio for both neutrino and antineutrino 
running. 

For the 8 GeV running we have used a likelihood performance that has 80% efficiency for signal 
in figure [7j Such a choice uses the advantages of running a low energy beam which will have lower 
backgrounds from both the neutral currents and the beam contamination. The advantages of the 
lower energy beam can be clearly seen in the large CP effects below 1 GeV. The appearance signal 
is seen to vary by more than a factor of 2 over the range 5cp = — 7r/2 — )• tt/2. It is clear that for 
the 8 GeV beam, the antineutrino event rate will be low. The principal reason for the low rate 
is the antineutrino charged current cross section at low energies. In particular, the quasi-elastic 
cross section for antineutrinos reduces much more rapidly at low energies than for neutrinos. The 
second reason is the lower production of tt~ mesons by lower energy protons. It should also be 
noted that the large CP effects below 1 GeV are only weakly affected by mass hierarchy. 

A simultaneous run of beams at 60 GeV and 8 GeV will provide excellent coverage of appearance 
spectra across the entire energy region. Moreover, the beam running could be managed in such 
a way that the two beams provide opposite polarity beams. For example, the 8 GeV beam could 
run in the neutrino mode and the 60 GeV beam could run in the antineutrino mode. In the case 
of inverted hierarchy, such a mode of running could be extremely beneficial. 

We have used the spectra shown in figures [6] and [7] to calculate the precision with which 
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FIG. 7: Event spectra using a 200 kTon water Cherenkov detector at Homestake for a neutrino beam 
from FNAL. Wc have assumed 3 MW of power at 8 GeV for 5 years of running for neutrinos (top plots) 
and antineutrino modes (bottom plots). The left plots are for normal mass hierarchy and right plots are 
for inverted mass hierarchy. We have put sin 20i3 = 0.1. The curve with error bars in each plot is for 
Sep = — 7r/2 and the red and blue curves are for Sep = 0, tt/2, respectively. The integral numbers of events 
are shown in the legends. The integral was computed up to the arrow in the figure. 



sm'^29is,Scp could be measured. These calculations were performed using the GLoBES software 
program [TU]. We have investigated the sensitivity with respect to moving the energy threshold 
From 100 MeV to 500 MeV and found the sensitivity to be only weakly dependent on this thresh- 
old. The precision is shown in figures [8] and [9j For this plot we have used a total of 5 yrs of 
running in 60 GeV and in 8 GeV beams (as remarked earlier this could be simultaneous). The 
60 GeV running is split in neutrino and antineutrino modes, but the 8 GeV running is in only 
the neutrino mode. For the combined plot, we have also included the data from the pre-project-X 
period for LBNE which would span a total of 10 years with 700 kW of running at 120 GeV |11) . 
The precision is dominated by the 60 GeV running after Project-X, nevertheless the pre-Project-X 
running certainly helps. 

We will remark on a few features of this measurement: 



• The two beams cover very different energy regions and have independent sensitivity. The 60 
GeV data will be affected by large matter effects and will result in the resolution of the mass 
hierarchy with very high (> 15cr) confidence. Once the mass hierarchy is known, the same 
data can be used for additional constraints on the matter potential or other contributions 
to the oscillation as described below. 

• The 60 GeV data provides high precision on sin^ 2^13 (~ few percent). The 8 GeV data with 
neutrino only running will have high precision to dcp, but it will be correlated to ^13. When 
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the two sets of data are combined, a measurement of 5cp with an error of ±10° (±5°) at 
6cp = 7r/2 {6cp = 0) is possible. 

• It is certainly possible to run the 60 GeV beam for much longer running time and obtain 
approximately the same precision on {613, Scp) as the combined 8 GeV and 60 GeV result. 
However, this will not provide the same independent constraint on the parameters that could 
lead to discovery of new phases or interactions. 

• The figures |8] and [9] were computed for normal mass hierarchy. For the inverted mass hi- 
erachy the sensitivity is the same because the statistical merit of the combined neutrino 
and antineutrino data set for the 60 GeV beam remains approximately the same. For the 
inverted hierarchy there will be higher signal in the antineutrino mode and lower in the 
neutrino mode. However, it should be noted that for the inverted hierarchy the 8 GeV beam 
should still be run in the neutrino mode. 

• The two independent data sets form independent constraints on the neutrino oscillation 
parameters. If there is any additional potential difference between the two mass eigenstates 
due to new physics it will become evident as a shift in the measured parameters from these 
two datasets at different energies. Such a shift could appear even if there is no explicit CP 
violation in the 3-generation picture as represented by the parameter Scp- This measurement 
precision could be used as a model independent parameter to test for new physics. 

• All calculations of sensitivity so far have assumed that the true value for sin 2^23 = 1-0 or 
^23 = vr/4. However, 623 is actually not very well known. The measurement ranges from 
35.7° to 53.2° at 3a. In particular, the value could be above or below 45°. The deviation of 
^23 from maximum mixing is of high interest to GUT theorist who can predict this particular 
mixing angle through various models. The deviation of ^23 from 45° might be indicative of 
the relationship between quark and lepton mixing. 

Our current calculations include a precise determination of sin'^ 2^23 using the disappearance 
data from LBNE. This internal constraint is included in the sensitivities published for LBNE. 
However, the disappearance data is not sensitive to the octant of 023- The appearance data 
has excellent sensitivity to the angle ^23 as shown in equation [T] The first term of the 



probability depends on sin^ ^23- This can be seen in figure [lOj The oscillation probability at 
high energies using the 60 GeV beam is seen to have a large dependence on 623, but the 8 
GeV spectra do not display the same sensitivity. A joint fit using 8 GeV and 60 GeV data 
is expected to resolve the octant of ^23- We will display this calculation in an update to this 
note. 



Conclusion 

We have examined a method to obtain multiple constraints in the neutrino oscillation sector 
using Project-X with high statistics, low energy neutrinos over long distances. This is motivated 
by the recent indications of a modest value for ^13 ~ 9°. For oscillation physics at low energies 
the charged current cross section is dominated by quasielastic scattering. For the quasielastic final 
state any detector capable of measuring single lepton final states is adequate, however as a result 
of the work on LBNE the water Cherenkov detector has been shown to be cabable of achieving the 
needed target mass to obtain the statistical precision. Furthermore, LBNE related simulation and 
reconstruction work has resulted in good understanding of the performance of this detector in the 
beams that we have described above. And therefore we have used a 200 kTon detector to obtain 
event spectra and sensitivity in the {0i3,6cp) space. 

Another strategy to obtain events at high values of L/E is by utilizing a much longer baseline 
such as 2500 to 2700 km from FNAL to California or Oregon. We suggest a renewed examination 
of such a strategy and optimization of the beam spectrum for much longer baselines. For example. 
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FIG. 8: 2 sigma and 1 sigma confidence level measurements of sin 29i3,6cp with different running con- 
ditions. Left-most plot is for running 60 GeV beam for 2.5 yrs each in neutrino and antineutrino modes. 
Middle plot is for running the 8 GeV beam only in the neutrino mode for 5 years. The right-most plot is the 
combination of both running conditions combined with preproject-X running of 10 yrs with 120 GeV and 
700 kW beam. For this calculation we have assumed sin^ 20i3 — 0.1, Sep = 7i'/2 and normal mass ordering. 
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FIG. 9; 2 sigma and 1 sigma confidence level measurements of sin 29i3,Scp with different running con- 
ditions. Left-most plot is for running 60 GeV beam for 2.5 yrs each in neutrino and antineutrino modes. 
Middle plot is for running the 8 GeV beam only in the neutrino mode for 5 years. The right-most plot is the 
combination of both running conditions combined with preproject-X running of 10 yrs with 120 GeV and 
700 kW beam. For this calculation we have assumed sin^ 26i3 = 0.1, 5cp = and normal mass ordering. 

running with different tunes of the FNAL neutrino could allow good coverage of L/E over much 
longer baselines. 

Project-X will allow great flexibility in proton beam energy and power that can be delivered 
for neutrino beams. In particular, we find the possibility of an upgrade to the pulsed LINAC for 
Project-X to deliver high power at 8 GeV very interesting. With such an upgrade high power 8 
GeV running can be simultaneous with 60 GeV running because only 266 kW of the 8 GeV beam 
will be fed to the Main Injector to make 2 MW of 60 GeV. The rest of the almost 4 MW of 8 GeV 
could be used separately for low energy neutrino production. A combination of 8 GeV and 60 GeV 
running will allow coverage of multiple oscillation nodes with high statistics. 

The combined data set from 8 GeV and 60 GeV running will allow multiple independent con- 
straints on the mixing parameters. The shift in the 9i3,6cp solutions from these data sets could 
be used as a model independent parameter to test for new physics. If the data sets match, very 
high precision ranging from ±5° to ±10° on 6cp is possible with a few years of running. 
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FIG. 10: Effect of 6*23 variation on tlie appearance spectra from 8 GcV running (left) and 60 GeV running 
(right). The integral numbers of events are shown in the legends. The integral was computed upto the arrow 
in the figure. 
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